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Abstract 

The electronic structure of the infinite layer compound CaCu02 has been calculated with the B3LYP hybrid density func- 
tional. The mixing of the Hartree-Fock exchange in the exchange-correlation energy separated the bands crossing Fermi energy 
to form an antiferromagnetic insulating ground state of charge transfer type. The complete elimination of the self-interaction 
through the exact exchange and the optimized gradient-corrected correlation energy significantly improved theoretical results. 
The theoretical energy gap and magnetic moment are in excellent agreement with the experiments. The ratio of intralayer to 
interlayer magnetic coupling constants and lattice parameters are also in good accordance with the experiments. Some char- 
acteristics of the electronic structure of insulating Sr2Cu02Cl2 from angle-resolved photoemission experiments are observed in 
the B3LYP band structure for CaCu02. 

PACS numbers: 71.15.Mb; 71.27.+a; 74.25.Jb 



As a parent compound of high temperature supercon- 
ductors(HTSC), the infinite layer material CaCu02 0] 
has a simple structure and a high transition tempera- 
ture when doped to optimum hole density]^, [sj. Similar 
to La2Cu04 and YBa2Cu30g, CaCu02 has a insulat- 
ing three-dimensional antiferromagnetic(AFM) ground 
state[Q, with energy gap A — 1.5eVQ and magnetic 
moment \x = 0.51 /isQ. Although the ratio of interlayer 
to intralayer magnetic coupling is one order of magni- 
tude higher than that in La2Cu04 and YBa2Cu306, the 
ratio is still very low. The nuclear magnetic resonance 
experiment on Cu in Cao.85Sro.i5Cu04 shows that the 
interlayer coupling is about two orders of magnitude less 
than the intralayer one|^. 

The electronic structure of CaCu02 has been investi- 
gated by several theoretical methods, such as LAPW0, 
LMTO-ASAg, FLMTOg. These LSDA-based meth- 
ods failed to give the correct ground state, as in the case 
of other HTSCs [|l^. Later, it was found that the fail- 
ure was due to the self-interaction inherent in the LSDA 
density functional, which tends to delocalize the elec- 
trons. The self-interaction correction(SIC) method gave 
correct AFM insulating ground state with /i = 0.58 
and A = 0.84 eV[0. The LSDA+C/ method also gener- 
ated correct ground state, with A = 2.1 eV and /i = 0.66 
\XB for U = 7.5 eV]l|l and A = 1.96 eV and [i = 0.71 [Ib 
for J7 = 5 eV|l^]. Although the two approaches recovered 
the correct ground state, their quantitative comparisons 
with the experiments were not satisfactory. 

In this paper we study the electronic structures of 
CaCu02 with the so-called B3LYP hybrid density func- 
tional. The hybrid functional was very successful in 
the thermochemistry of atoms and molecules . Later, 
the hybrid functional was applied to some periodic 
systems! 15, Qa]. The argument for mixing the Hartree- 



Exc is based on the adiabatic connection formula fl^. 
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where A is the electron-electron interaction parameter, 
with A = 1 corresponding to the real interaction and 
A = a system of noninteracting electrons. Uxc is the 
potential energy of the exchange-correlation at any A. In 
the case of A = C/^^ is the HF exchange energy of the 
corresponding noninteracting system. The weights for 
the gradient-corrected correlation energy, local exchange 
energy and the exact HF exchange terms were deter- 
mined by a linear least-square fitting of the thermochem- 
ical properties of some atoms and molecules. The atom 
with highest atomic number used in the fitting is CI. No 
atoms with d or higher shells were used. 20% mixing of 
the exact HF exchange energy in the Exc gives theoret- 
ical results in good agreement with experiments. In the 
B3LYP scheme the Perde w- Wang gradient-corrected 
correlation energy is replaced by Lee- Yang-Parr correla- 
tion energy jisf. 

The final exchange-correlation energy functional reads 
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Fock(HF) exchange in the exchange-correlation energy 



in which the local spin density functional of Vosko, Wilk 
and Nusair^ is used for Ej^^^^ and Ejj^^"^. Ef"^ 
is the exact nonlocal HF exchange energy. /S.E^^^ and 
AE^^^ are the Becke'spO| and Lee- Yang-Parr's gra- 
dient corrections for the local exchange and correla- 
tion energies, respectively. The optimum values for the 
parameters oq, ax and ac are 0.20, 0.72, and 0.81, 
respectively 0. 

In this paper the calculations are carried out with 
CRYSTAL package pl[. The basis vectors for expanding 
the Kohn-Sham orbitals are Bloch functions composed 



of localized contracted basis sets^^l- All-electron basis 
sets for Ca, Cu and O ions are of the form of 86-511G, 
86-411(41d)G and 8-411G, respectively. 

In the calculations 75 points in the irreducible part of 
the first Brillouin zone(FBZ) were used. We adopt the 7, 
7, 7, 7 and 14 as the integral tolerances to obtain high pre- 
cision in monoelectronic and bielectronic integrals. Also 
a strict criteria for convergence, i.e. when the differ- 
ence between two consecutive root-mean-squared values 
of the density matrix elements is less than 10~^ the con- 
vergence is assumed, is used to make the magnetic mo- 
ments on Cu sites converge. Different magnetic states, 
i.e. metallic(without spin polarization), ferromagnetic, 
two-dimensional(2D) and three-dimensional(3D) antifer- 
romagnetic structures, were checked to find the ground 
state magnetic configuration. For comparisons, the 
unrestricted Hartree-Fock(UHF) and general gradient- 
corrected LSDA(LSDA-I-GGA) calculations were also 
carried out with the same basis sets. We also performed 
structural optimization for the 3D AFM supercell. The 
optimization is implemented with Domin algorithm and 
the P4/mmm symmetry is preserved. 

Our results show that the 3D AFM state has the lowest 
energy, i.e. the B3LYP hybrid functional could predict 
correct ground state magnetic configuration. The energy 
difference per unit cell between the ferromagnetic and the 
2D AFM(antiferromagnetic order in the Cu02 plane and 
ferromagnetic order in the Z direction) ground state is 
0.285 eV. The 2D and 3D AFM states has nearly the same 
energies, with the energy of 3D AFM state is only 0.00164 
eV lower than that of 2D AFM state. The energies of 
different magnetic configurations could be described by 
the Ising model, 

E ^ ^ JijSiSj. (3) 

The intralayer coupling constant J|| and interlayer cou- 
pling constant J± could be obtained from the energy 
differences of the three different magnetic states. One 
can easily get J\\ — [EpM — E2dafm]/4: and J± = 
[E2DAFM - E3dafm]/2, in which Efm, E2DAFM and 
E3DAFM are the energies of FM, 2D AFM and 3D 
AFM states, respectively. The theoretical results are 
J|| = 825.5K and J± = 9.49K. The ratio J±/J|| is 1.1%, 
which is in good agreement with the experimental result 
that the interlayer magnetic coupling is about 2 orders 
of magnitude weaker than the intralayer couplingQ. 

The optimized lattice parameters for the supercell are 
a=5.563A and c=6.536A. The results are about 2% larger 
than the experimental values Q. The precision is compa- 
rable to LDA computations for conventional materials. 
In addition. Table || shows that the B3LYP theoretical 
magnetic moments and energy gap obtained with exper- 
imental and theoretical lattice parameters are in excel- 
lent agreement with the experiments ^. As mentioned 
above the optimum parameters were obtained with atoms 
and molecules without d-electrons, so the excellent agree- 
ment between the theory and experiments was unex- 



TABLE I: The energy gaps A(in eV) and magnetic moments 
//(in /ifl) from different theoretical schemes are compared with 
experiments for CaCu02 and La2Cu04. The results with su- 
perscripts a and b are obtained using experimental and opti- 
mized lattice parameters, respectively. 
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pected because of the strong localized character of d 
orbitals. But from the study of pseudopotentials one 
knows that the 2p orbitals in the first-row atoms also 
pose problems, these orbitals have also localized charac- 
ter because of the lack of orthogonal repulsions. So, if 
there are enough variational freedoms in the B3LYP hy- 
brid functional to correctly describe atoms with 2p shell 
then it could be applied to systems with 3d electrons. 

The large energy gap of the UHF results from the lack 
of correlation, i.e. the interaction is not screened. As can 
be seen in Table | that the UHF gaps for both CaCu02 
and La2Cu04 are about one order of magnitude larger 
than the experimental results. It is also much larger than 
the usual cases, which means that the correlation effects 
are more important in these materials. 

The projected densities of states(DOSs) are shown in 
Fig. 1 and the band structures along the high symmetry 
lines are shown in Fig. 2(a). These results are obtained 
with experimental structural parameters. The results 
show that CaCu02 is a charge transfer insulator. This 
is similar to the case of transition metal monoxides [p^. 
As shown in Fig. 1 the main components at the top 
of the valence bands and the bottom of the conduction 
bands are Cu 3da;2_j^2 and O 2p^(,y). Though the spectral 
weight at the bottom of the conduction bands is mostly 
of Cu 3(1^2 _y2 nature, there is a significant ingredient of 
O 2p2.(y) at the top of the valence bands. Because of 
the lack of apical oxygens in CaCu02 the bonding along 
the Z direction is weak, which is different from La2Cu04, 
where significant Cu 3(1^2 / O 2pz components are existent 
at the top the valence bands [|l6| . The result suggests that 
the Cu 3dz2/0 2p2 components in La2Cu04 may not be 
essential to the hight temperature superconductivity of 
the system. 

One may also notice that some features in Fig. 2(a) 
are in good agreement with the angle-resolved pho- 
toemission(ARPES) experiments on the insulating 
Sr2Cu02Cl2Q. Sr2Cu02Cl2 has a structure similar to 
La2Cu04, with the apical oxygens in La2Cu04 replaced 
by Clp^. The X and M points in Fig. 2(a) correspond 
to the (7r/2, 7r/2) and (tt, 0) points in the FBZ of the 
Cu02 plane, respectively. The (tt, tt) point in the ARPES 
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experiments is equivalent to F point in Fig. 2(a). The 
ARPES experiments show that the valence band top is 
at (7r/2, 7r/2) and there is a nearly isotropic dispersion 
aromid this point. The dispersion from (tt, 0) to (0, 0) 
is flat and the dispersion from (tt, 0) to (0, tt) is similar 
to the one from (0, 0) to (tt, tt). All these character- 
istics are consistent with om" results, as can be seen in 
Fig. 2(a). Due to the folding-up of the FEZ in the case 
of the AFM supercell, there are two nearly degenerate 
bands around X points. Although the dispersion along 
the F-X-M direction in Fig. 2(a) is stronger than experi- 
mentally observed, a better agreement with experiments 
may be expected by applying the method to Sr2Cu02Cl2. 

To see the effect of the mixing of the exact exchange 
the LSDA-t-GGA band structure is shown in Fig. 2(b) for 
comparison. The same supercell and same functionals for 
GGA corrections were used in the LSDA-f GGA and the 
B3LYP calculations. From Fig. 2(a) and Fig. 2(b) one 
can see that the mixing of the exact exchange results 
in no significant change of the dispersions and relative 
positions of the bands far from Fermi energy Ep. But 
the important effect of mixing exact exchange is that two 
bands, which cross Ep, are separated from the other two 
bands. The separated two bands have no overlap with 
the other conduction bands, which is different from the 
LSDA-hf/ results[|l|, In the LSDA-f [/ scheme there 
are also two separated bands composed of Cu 'idx2_y2 
and O 2]ix{y) antibonding states, but the two bands are 
overlapped with other conduction bands p^. 

Except the difference mentioned above between the 
B3LYP and the LSDA-I-C/ schemes, some essential char- 
acters are shared by the two approaches. The reduced en- 
ergy gap results from upward shift of part of Cu 3da;2_j,2 
and O 2p3;(j,) spectral weights, and this part of spec- 
tral weights have similar width in both schemes. Also 
the most important projected DOSs for Cu 3d2.2_j^2 and 
O 2p3,(y) have the same characters . Although the 
LSDA-|-[/ captures the essential physics in strongly cor- 
related systems its quantitative results are not as good as 
the B3LYP. The reason may be due to the simple mean- 
field treatment of the Hubbard U term in the LSDA-|-[/. 
Recently, Mazin et al. have shown that the LSDA-I-C/ 
failed for moderately-correlated metal, such as FeAl|p3|. 

In the LSDA-t-f/ scheme it is easy to see the reason 
why the bands crossing Ep were separated to form a 
gap. The additional orbital dependent potential splits 
a band into Hubbard subbands, with separation equal 
to approximately the screened on-site Coulomb repulsion 
U. The SIC approach restores the localized character 
of d orbitals, thus introducing strong on-site Coulomb 
repulsion. In the studies of the Mott insulators it was 
found that the simple UHF can give qualitatively correct 



results, whereas more advanced DFT methods failed p4[. 
The reason is that in the UHF there is no self-interaction. 

In our UHF calculation for CaCu02 the energy gap is 
about one order of magnitude larger than the experimen- 
tal result. It's expectable that small amount of mixing 
of the exact exchange would reduce the unphysical large 
gap to the experimental value while maintaining the cor- 
rect ground state magnetic configuration. But in the 
B3LYP approach it is not very clear why only 20% mix- 
ing of the exact exchange into the GGA corrected LSDA 
can greatly change the LSDA-f GGA results. To see the 
point one can rewrite Eq. (0) as 
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Instead of only 20% exact exchange one can view the 
B3LYP as incorporating 100% exact exchange, and treat 
the second and the third terms as additional contribu- 
tions to the correlation energy. So, with the exact ex- 
change, where the self-interaction is completely excluded, 
and additional variation freedoms for correlation energy 
one should expect that the method would give better re- 
sults than the LSDA-f GGA and the UHF approaches, in 
the latter case the correlation energy is totally discarded. 

One may view the B3LYP hybrid functional as a much 
improved energy functional towards the exact one. For a 
specific class of materials better results can be obtained 
by reoptimizing the parameters appearing in the corre- 
lation energy and by adding additional functionals. It is 
also interesting to test if the B3LYP functional could be 
successfully applied to moderately-correlated and doped 
strongly correlated systems. Muscat et al. have applied 
the B3LYP to materials of different types of bonding [p^. 
The theoretical energy gaps are in good agreement with 
experiments, and the band structure for Si is also in good 
accordance with experiments and other theoretical meth- 
ods, such as quantum Monte Carlo. For the Mott in- 
sulator NiO the B3LYP also generated much better en- 
ergy gap than the LSDA-f C/[|| and the SIC-LSD|2|, |6). 
So, one can expect more successful applications of the 
method to other materials. 

To conclude, the B3LYP scheme has been applied 
to CaCu02. The mixing of the exact exchange elimi- 
nated the self-interaction and caused the band separation 
about the Fermi energy compared with the LSDA-f GGA 
method. The magnetic moment and energy gap are in ex- 
cellent agreement with experiments. Some of the struc- 
tural and magnetic properties are also in good agree- 
ment with experiments. The ARPES characteristics on 
Sr2Cu02Cl2 are observed in our calculations on CaCu02. 
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Fig. 1 The projected densities of states of Cu 3d, O 
2p and Ca 3sp partial waves in a V2 x V2 x 2 supercell 
of CaCu02. All the DOSs are for spin-up electrons and 
the magnetic moment of Cu is -0.51/iB. 



Fig. 2 (a) . The B3LYP energy bands of a \/2 x \/2 x 2 
supercell of CaCuOs. (b). The LSDA-f GGA band struc- 
ture of the same supercell. The coordinates of the high 
symmetry points are F = (0,0,0)/2, M = (l,l,0)/2, 
Z = (0,0,l)/2, R = (0, l,l)/2 and A = (l,l,l)/2. The 
bands are for spin-up electrons. 
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Figure 2(a) 
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Figure 2(b) 
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